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Abstract. We investigate the connection of lattice calculations of moments of isovector parton distributions 
to the physical regime through extrapolations in the quark mass. We consider the one pion loop renor- 
malisation of the nucleon matrix elements of the corresponding operators and thereby develop formulae 
with which to extrapolate the moments of the unpolarised, helicity and transversity distributions. These 
formulae are consistent with chiral perturbation theory in the chiral limit and incorporate the correct heavy 
quark limits. In the polarised cases, the inclusion of intermediate states involving the A isobar is found to 
be very important. The results of our extrapolations are in general agreement with the phenomenological 
values of these moments where they are known, and for the first time we perform an extrapolation of the 
low moments of the isovector transversity distribution which is consistent with chiral symmetry. 

PACS. 12.38.Gc - 11.30.Rd 



1 Introduction 



Until recently, state of the art studies Oj] of the moments of 
parton distributions within lattice QCD have led to large 
. discrepancies with experimental data Lattice calcula- 
tions are performed at quark masses much greater than 
those of the physical light quarks. Consequently, an ex- 
trapolation to the physical mass regime must be made 
in order to compare with experiment. A naive linear ex- 
trapolation results in values for the first three non-trivial 
moments of the unpolarised u — d distribution that are 
50% above the phenomenological values For the po- 
larised distributions the results are no better; (the th 
moment of Au — Ad) is significantly underestimated by a 
linear extrapolation, while the lack of accuracy in the data 
(both experimental and lattice) on higher moments pre- 
cludes definitive statements. Such a large disagreement in 
such basic hadronic observables casts doubt on the current 
reliability of the lattice approach to hadronic physics. 

We present an improved extrapolation scheme ||,|],[|, 
|(| that for the first time resolves much of this discrepancy. 
Using constraints from chiral symmetry and the heavy 
quark limit, we develop a formalism for the extrapolation 
of the moments of the isovector, unpolarised, helicity and 
transveristy distributions. These are related to the for- 
ward nucleon matrix elements of various twist-2 operators 
which are calculated on the lattice through the operator 
product expansion. In particular we find that contribu- 
tions from intermediate states involving the A isobar are 
large, and their effects cannot be ignored in any quanti- 
tative analysis. With their inclusion, we are able to make 



reliable, almost model-independent^ extrapolations of the 
moments of parton distributions and make predictions for 
the low spin transversity moments which can be tested by 
future measurements. 

2 Pion-Loop Renormalisation of Moments of 
Parton Distribution Functions 

General constraints from the approximate chiral symme- 
try of QCD lead to the appearance of non-analytic terms 
in the quark mass expansion of many hadronic quantities. 
In particular, the moments of quark distributions behave 
as §| 

(x n ) q ,Aq,Sq ~ TO 2 log TO 2 , (1) 

where the Gell-Mann-Oakes-Renner relation, m q ~ to 2 , 
has been used to express the quark mass in terms of the 
pion mass. This behaviour arises from the infrared prop- 
erties of the one-pion loop renormalisation of the matrix 
elements of the corresponding twist-2 operators O q< Aq,Sq 
(see Ref. H for their definitions and further details). That 
is, 

(N\0^ n \N) dressed - ^(jV|0f-^|iV) baxe) (2) 

where Zi is the nucleon wavefunction renormalisation and 
Zi, i = q, Aq, 5q, are the operator renormalisations arising 
from irN, ttA, and ir(NA transition) intermediate states. 

1 Results are independent of the shape of the pion-nucleon 
form factor to 1-2% §]. 
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Fig. 1. Lowest three nontrivial moments of the unpolarised 
isovector parton distribution. Data are taken from various lat- 
tice simulations (see Ref. || for details). The linear extrap- 
olation (dashed) significantly overestimates the experimental 
results (stars), whilst the extrapolation using Eq. (0) is in rea- 
sonable agreement. 



Fig. 2. Lowest three nontrivial moments of the helicity isovec- 
tor parton distribution. Curves shown are (dotted line) linear 
extrapolation, (dashed line) extrapolation using the polarised 
analog of Eq. (ph but ignoring A contributions, and (solid lines) 
extrapolations for varying values of the tvNA coupling. (Data 
as for Fig. |l|) 



In the case of the unpolarised moments, simple one- 
loop calculations of the renormalisation using a variety 
of form factors (cutoff, monopole, dipole) for the pion- 
nuclcon coupling all lead to a simple extrapolation formula 



1 + c L na?t4 log ■ 



with three free parameters, a r , 



(3) 



/« (clna is fixed by 



chiral perturbation theory). Similar expressions are found 
for the helicity and transversity moments in Ref. [pi . The 

§ameter \i describes the physical scale of the pion source 
and ideally would be constrained by lattice data. How- 
ever, until sufficiently accurate lattice data at low quark 
masses are available, alternative methods (such as the one 
used here) must be used to fix /i. 

As shown by the solid curves in Fig. |l| this formula 
provides an excellent description of the lattice data, the 
experimental values of the moments and (with suitable 
modification - see Ref. ^]) the values known in the heavy 
quark limit. However, when a similar formula is applied 
to the moments of the helicity distributions, the extrapo- 



lation is considerably worse than a naive linear extrapo- 
lation, as shown by the dashed lines in Fig. [|. 

When we turn to the moments of polarised parton dis- 
tributions, there is considerable evidence from phenomeno- 
logical models that suggests the A resonance will play an 
important role. Although the A contributions formally en- 
ter at higher order in m w , the coefficients of these next- 
to-leading non-analytic terms are large, and they cannot 
be ignored in any quantitative analysis. This is clearly 
demonstrated by Fig. ^ where the pion mass dependence 
of the one- loop renormalisation of helicity matrix elements 
is shown calculated with a dipole form factor with varying 
values of the ttNA coupling, g^NA- Full details are given 
in Ref. ||. For these helicity (and transversity) operator 
matrix elements the difference between g-nN a/ g^N n = 
(no A) and g-nN a/ 9-kNN = 2 (phenomenological value) is 
up to 15% — much greater than in the unpolarised case 
where the overall effect is less than 2% over the entire 
mass range studied here. 

In order to incorporate these effects into the extrap- 
olations of lattice data, we calculate the required renor- 
malisations for the phenomenologically preferred dipole 
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Fig. 3. One-loop renormalisation of the nucleon matrix el- 
ements of the polarised, isovector, twist-2 operators using a 
dipole form factor (A = 0.8 GeV) and varying values of the 
coupling ratio g-jrNA/gwNN- The shaded region corresponds to 
variation of the Weinberg- Tomozawa coupling - see Ref. Q|. 

Moment Experimental Extrapolated 

0.145(4) 0.17(3) 

0.054(1) 0.05(2) 

(x 3 )„_ d 0.022(1) 0.02(1) 

(l)Au-Ad 1.267(4) 1.12(8) 

(x) A u-Ad 0.210(25) 0.27(3) 

(x 2 )Au-Ad 0.070(11) 0.14(4) 

(l)i„-«d ? 1.22(8) 

(x) 5u - 5d ? 0.5(1) 

Table 1. Values of the unpolarized, helicity and transversity 
moments, extrapolated to the physical pion mass using Eq. (^). 
For comparison, experimental values of the moments, where 
known M, are also listed. 



parameter, A = 0.8 GeV, and for varying values of the 
coupling ratio g^NA/g^NN- We then determine ji by fit- 
ting to the calculated renormalisations using Eq. (0) (with 
Clna calculated in the AM = M A - M N -> limit). With 
/X thus fixed, we then use the lattice data to determine the 
fit parameters a n and b n . The resulting curves are then 
shown in each panel of Fig. || for g^NA/g^NN = (no A), 
v/72725, 1.85 and 2.0 

Table [j] shows the resulting extrapolated values along 
with uncertainties resulting from the irNA and Weinberg- 
Tomozawa couplings, statistical and (estimated) system- 
atic errors. Also listed are the experimental moments where 
known (2| . It is evident that the agreement between the ex- 
perimental and extrapolated results is very good for the 
unpolarised moments, but significantly worse in the he- 
licity sector. The extrapolated value of gA is 10% below 
its experimental value, with 8% errors. However, there is 
some evidence that lattice simulations of this quantity are 
particularly susceptible to finite volume effects 0|. For 
the higher helicity moments, statistical errors on the lat- 
tice data are not yet sufficiently small to make definitive 
statements. Since the renormalisation of the transversity 



matrix elements is almost identical to that of the helicity 
matrix elements, we are also able to make an estimate of 
the values of these moments which can be compared with 
future experimental determinations. 



3 Conclusion 

We have investigated the extrapolation of lattice data on 
the low moments of the isovector, unpolarised, helicity and 
transversity parton distribution functions. In the case of 
the polarised moments, the A isobar, even though not 
leading non-analytic, was found to play a significant role 
and its inclusion in the operator renormalisation and ex- 
trapolation procedure is vital. With this accounted for, 
we have the surprising result that the effect of the non- 
analytic behavior is strongly suppressed for the polarised 
moments, and a naive linear extrapolation of the moments 
provides quite a good approximation to the more accurate 
form. 
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